New Dissipation Relaxation Phenomenon in Oscillating Solid 4 He 
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We describe the first observations on the time-dependent dissipation when the drive level of a 
torsional oscillator containing solid 4 He is abruptly changed. The relaxation of dissipation in solid 
4 He shows rich dynamical behavior including exponential and logarithmic time-dependent decays, 
hysteresis, and memory effects. 

PACS numbers: 67.80.-s, 61.72.Hh, 62.40.+i 



The recent discovery [l|, Q and its confirmations 
[EH, [EH] °f significant low temperature reductions in the 
rotational moment of inertia of solid 4 He below 300 mK 
have generated enormous excitement as possible evidence 
for the long-sought supersolid behavior [jj. The physi- 
cal origin of the apparent non-classical rotational inertia 
(NCR!) extracted from the temperature dependent shift 
in the torsional oscillator frequency, however, remains un- 
settled Q. The majority of the torsional oscillator exper- 
iments showing NCRI has been devoted to measuring the 
oscillator frequency under a variety of conditions in such 
as crystal quality [3J, sample size and geometry, and 3 He 
impurity concentration Relatively little experimental 
information is available on the dissipation accompanying 
the torsional oscillation. If the observed reductions in 
rotational inertia are related, for example, to vortex liq - 
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or defect motion [l4[ in solid 4 He, then dynamical time- 
dependent effects might occur in the dissipation in the 
oscillating solid 4 He samples. We carried out systematic 
observations on the time-dependent dissipation accompa- 
nying sudden changes in the torsional oscillator displace- 
ment amplitude. We have discovered that the dissipa- 
tion in solid 4 He shows rich dynamical behavior including 
exponential and logarithmic time-dependent relaxation, 
hysteresis and memory effects. 

The torsional oscillator utilized was the compound os- 
cillator identical to that described in Q. Briefly, our 
oscillator contained an upper "dummy bob" hung by an 
upper rod and connected to a lower rod attached to the 
sample chamber containing a cylindrical (10 mm diame- 
ter and 8mm height) solid 4 He . The oscillator resonated 
at "in- phase" and "out-phase" modes at frequencies, /i 
= 496 Hz and = 1173 Hz, respectively, whose ampli- 
tudes and phases were continuously tracked. The tem- 
perature was measured by a calibrated 3 He melting pres- 
sure thermometer mounted onto the vibration isolation 
copper block to which our oscillator was attached. All of 
our solid samples were grown from commercially avail- 
able 4 Hc gas with nominal 0.3 ppm 3 He impurity by the 
usual blocked capillary method. 

The procedure for measuring the time dependent be- 
havior of our torsional oscillator is illustrated in Fig. [1] 
The f% mode resonance was continuously tracked by a 
lock-in amplifier maintaining a constant phase between 
the drive and the response. The shifts in resonance fre- 
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FIG. 1: Time dependence of, ac drive amplitude (a), fre- 
quency shift (b), and rim velocity ((c) and (d)) at T = 9.95 
mK. The initial v r (~160 /im/s) corresponds to displacement 
amplitude of 22 nm. The drive (torque) level is tx VdcV ac , 
where the dc bias voltage Vdc is kept constant. The ring 
down portion of v r is not shown in (d). The drive level was 
increased back to the original level at t = 5xl0 4 s. The fre- 
quency shift does not return to the initial value in accordance 
with our observation of hysteresis After the increase in 
drive, there is no decay similar to that after the first decrease 
in drive. The ring down time of the empty oscillator itself 
ranges from 120 to 200 s. 



quency and the response amplitude were monitored con- 
tinuously throughout the procedure. The response am- 
plitude was converted Q to sample rim velocity ampli- 
tude, v r . The cell temperature was initially raised to 
90 mK, where the NCRI was found to be reversible as a 
function of drivel level [|| . The oscillator was then driven 
at a given level and allowed to achieve a steady state. 
The temperature was next lowered within about 10 3 s 



to 9.95 mK at the time defined as t 
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s and subsequently regulated within 10 /iK. The drive 
level was reduced to half of the initial value at t = 
(see panel (a)). Following the reduction in drive, the rim 
velocity (panel (d)) unexpectedly "undershot" and even- 
tually recovered in a complex relaxation process. The 
usual ring down of the oscillator resonance transpired 
during the time interval prior to the recovery from the 
undershooting. When the drive level was restored to the 
original value at t = 4.4 x 10 4 s, v r increased more than 
twice but the relaxation behavior was absent. When 
the process was repeated to decrease the drive level by 
half at t = 5.4 x 10 4 s, the initial undershoot behavior 
was now absent but "remembered" its previous state at 
the same drive level. Thus, the unexpected relaxation 
phenomenon occurs only in the first decrease from one 
high drive level to a lower one. Majority of the relax- 
ation phenomena described here was studied with the f% 
mode. The fi mode was studied at 9.95 and 19.65 mK 
and qualitatively similar relaxation phenomena were ob- 
served. When the solid 4 He was replaced with superfluid 
and the similar procedure was followed, the undershoot 
and subsequent recovery were both absent at all temper- 
atures studied (see Fig. [3]). This provided clear evidence 
that the observed unusual time dependence in the oscil- 
lator response originates in the solid 4 He sample. 

During the observed decay process, the oscillator fre- 
quency also showed time dependence (see Fig. QJb)). 
When v r was initially decreased, the oscillator frequency 
(therefore, the NCR! fraction) increased. This is the 
" critical velocity" effect as seen previously [U, [|| . There 
was an associated relaxation in frequency shift which ap- 
peared with similar dynamical behavior as v r . Owing to 
the much smaller effect on the frequency relaxation, we 
focused on the much larger effects of dissipation dynam- 
ics. Observations on the frequency relaxation resulting 
from temperature changes were reported recently (TEj . 

To see how the amplitude relaxation phenomena as 
shown in Fig. [1] depended on temperature, the same pro- 
tocol was followed except the final temperature. Some 
of such measurements are shown in Fig. [5] Similar un- 
dershoot in rim velocity was seen at temperatures below 
50 mK just after the decrease in drive level. Following 
the relatively rapid recovery, there was a much slower re- 
laxation at t > 10 4 s. The long time tail behavior after 
the decrease of drive level was complicated. In the long 
time limit, the rim velocity decreased in the temperature 
range, 20 mK > T > 30 mK, but it increased outside of 
this range. 

There have been suggestions that the supersolid state 
might involve glassy behavior in some way. The decay 
of magnetization of spin glass systems [l6| has been ob- 
served under certain conditions to decay with stretched 
exponential time dependence. Fitting our data to such 
time dependence was attempted but inadequate. We find 
that a sum of purely exponential and logarithmic time 
dependencies can fit our data satisfactorily: 



where vq, B, C and r are fitting parameters. The short 
and long time recovery behaviors are described by the ex- 
ponential and the logarithmic time dependence, respec- 
tively. The parameters vq and B are first set by fitting 
the data in the long time limit (t > lxlO 4 s), where the 
exponential term is negligible. The "origin" of time is 
taken as that when the final temperature is stabilized 
(t = —to) to account for the decay during the waiting 
time interval to. Once the logarithmic decay parameters 
are set, the remainder is fitted with the exponential de- 
pendence in Eq. (fT]). The fits (lines) with parameters 
shown in Fig. [5] appear adequate. The exponential de- 
cay amplitude C decreases as the temperature is raised 
and vanishes in the fit at 66.6 mK in Fig. 
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v Q + B ln(i + t ) - C exp(-t/r) 
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FIG. 2: Time dependence of the sample rim velocity (v r ) at 
different temperatures. The same protocol as in Fig. Q]was 
followed. The temperatures were 9.95, 29.78, 47.85 and 66.60 
mK. The lines are fits to Eq. (fTJ discussed in the text. 

To examine the dependence of the relaxation behavior 
on the initial drive level, observations were made at 29.45 
mK in the same protocol as above keeping the same fi- 
nal v r at ~ 22 /zm/s but varying the initial v r between 
180 and 100 /im/s. As shown in Fig. [3J the undershoot 
magnitude increased as the difference between the initial 
and final v r was increased. A small dependence on the 
difference was also present (see below) in the long time 
tail behavior. When the cell was filled with superfluid 
4 He, there was no relaxation other than the oscillator 
ring down. The inset shows the dissipation in the sample 
(subsequent to the end of the undershoot) computed from 
the observed amplitude A (Aq) in the presence of solid 
(superfluid) 4 He sample via 6(1/ Q) = (1/Q) - (1/Qo), 
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FIG. 3: (color online) Time dependence of sample rim velocity 
for different initial drive levels at 29.45 mK. The initial drive 
level was set such that v r were 104 (circles), 143 (pluses) 
and 180 (triangles) /im/s and the same final v r « 22 /xm/s. 
When the cell was filled with superfluid, there was only the 
oscillator ring down with no anomalous undershoot response 
(dots). The inset shows the dissipation due to the presence 
of solid 4 He sample derived (see text)from the data for each 
initial v r in the main panel. 

where Q = aA and Qq = aAo (a is the measured calibra- 
tion constant). The inset demonstrates the unexpected, 
large transient excess dissipation observed during the re- 
covery process just after the change in drive level. 

Temperature dependence of the fitted time constant r 
is displayed in Fig. 21 which includes those from the mea- 
surements with the same protocol as in Fig. [1] but with 
data taken only up to about t = 2 x 10 4 s. As the tempera- 
ture is increased, r decreases monotonically. The temper- 
ature dependence may be represented by r = Toe A ^ kBT \ 
t = 5.4xl0 2 s and A/fee = 14 mK. This characteristic 
temperature represents the onset of a newly uncovered 
dissipation dynamics in oscillating solid 4 He. 

The logarithmic variation in the relaxation in Eq. ([1]) 
is a common time dependence to various decay phenom- 
ena in superfluids [l7|, superconductors frictional 
processes [l9| and glasses [l6j]. The logarithmic relax- 
ation we observe has complex dependence on tempera- 
ture and v r . The sign of the fitted parameter B is neg- 
ative around 30 mK for the particular protocol followed 
for the data in Fig. [2] It is positive (negative) when the 
final v r is greater (less) than rj 80 /im/s, but this rim 
velocity boundary is dependent on the initial v r as well 
as temperature. The effects of different sample prepara- 
tion procedures and annealing on the logarithmic decay 
have not yet been studied. Judging from the hysteresis 
behavior of NCRI fraction being relatively insensitive to 
different samples in our cylindrical cell, we expect similar 
insensitivity in the relaxation phenomena described here. 



The behavior of the parameter B is complicated. It 
is known that the particular thermal history ([l5| and 
our studies) of a solid 4 He sample affects the oscillator 
frequency shift and its relaxation. For example, the di- 
rection of decay depends on the direction of temperature 
change. We do not yet have a complete thermal history 
characterization of the logarithmic decay of the dissipa- 
tion and it is possible that the observed sign reversal of 
B is related to some subtle difference in thermal history. 

Possible association has been suggested between the 
observed NCRI phenomena and the dislocation line net- 
work depending on the number, n^, of condensed 3 He 
impurities on the lines 0, [l4[ • We consider the following 
effect of dislocation on our observed dissipation relax- 
ation phenomena. Increasing the drive level might par- 
tially breakaway 3 He from the dislocation lines. Then, 
the distance (Li) between adjacent pinning sites at 3 He 
locations increases [20]. Just after the drive level is de- 
creased at t = 0, tt-3 is initially too small for the new lower 
drive level and there is an excess dissipation resulting in 
the undershooting of the measured amplitude. As more 
3 He recondenses onto the dislocation lines, the dissipa- 
tion might decrease as seen by the subsequent increase in 
amplitude. A difficulty with this scenario is that assum- 
ing a uniform dislocation line density A ~ 10 5 cm~ 2 (20j 
(average separation between lines ~ 17 /im) and taking 
the diffusion coefficient D of 3 Hc impurities (of concen- 
tration X3) in solid 4 He as DX3 w 3 x 10 -11 cm 2 /s at 
low temperatures [2l[ give a diffusion time less than 10 
s, short compared to r in Fig. 01 The observed r might 
instead be related to the time for 3 He to migrate along 
dislocation lines towards an equilibrium distribution. 

Another possible mechanism for the observed relax- 
ation dynamics is the effect on the number of disloca- 
tions induced by the change in oscillation amplitude. 
Hiki and Tsuruoka [22j showed that suddenly imposed 
changes in temperature AT could exert internal shear 
stresses and thereby change the number of dislocations. 
The measured relaxation in ultrasonic attenuation was 
interpreted as evidence for quantum mechanical tunnel- 
ing of dislocations [13] • The tunneling rate increased as 
the imposed change in temperature was increased. If the 
change in oscillator amplitude in our experiment were 
equivalent to (but likely to be much less than) AT = 
(1/10)T, the resulting tunneling time constant would 
be 10 5 s, much longer than the observed value. The 
quantum mechanical tunneling mechanism also leads to a 
temperature-independent relaxation time contrary to our 
finding. Thus, it is not likely that the relaxation time we 
observe is related to the dislocation tunneling. The ob- 
served long-term memory effects and logarithmic time de- 
pendence are also difficult to explain by the dislocation- 
related mechanisms. 

We discuss next our observations in terms of a vor- 
tex liquid model of supersolid state. In the model, 
NCRI is interpreted as the rotational susceptibility of the 
quantized vortex tangle and the dissipation arises from 
the interaction between the torsional oscillation and the 
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FIG. 4: Temperature dependence of fitted time constant r in 
Eq. [1] This relaxation time does not depend on the drive 
level. 

thermally fluctuating vorticics. Prior to the first change 
in the drive level, our procedure of initiating the oscil- 
lation at a high temperature is likely to have brought 
the sample solid 4 He to a quasi-equilibrium array of vor- 
tices for the drive level at the temperature of measure- 
ment. The observed logarithmic time dependence proba- 
bly arises from the thermally activated motion of vortices 
from one metastable state to another. This process con- 
tinues throughout the measurement while the sample is 
maintained at the constant temperature. Just after the 
drive level is decreased, the existing vortices are out of 
equilibrium with the new drive level. The observed un- 
dershoot may be a transient excess dissipation owing to 
the "wrong" number of vortices. The vortices now must 



adjust to the new drive level by, say, moving out of the 
sample. The characteristic time involved in the macro- 
scopic motion of the vortices is then the observed time 
constant r. The motion presumably involves processes 
occurring both within the sample and at the surface 
boundaries. The absence of overshoot and exponential 
decay in both the oscillator response and NCR! fraction 
implies that the drive reversal does not replenish vortices 
into the sample. This indicates the presence of surface 
barrier which inhibits creation of vortices in the system. 
When the drive level is decreased again, the number of 
vortices is already near that appropriate for the velocity 
and there is no undershoot nor accompanying exponen- 
tial relaxation. The time constant becomes shorter as 
the vortices begin to move more freely. This vortex liq- 
uid state appears to occur above about 60 mK. When the 
vortex liquid "freezes," the time constant would diverge. 
Our sample at the lowest temperature apparently did not 
attain a totally frozen vortex state. 

In conclusion, we studied the dynamics of dissipation 
imparted by cylindrical solid 4 He samples onto a torsional 
oscillator by suddenly decreasing the oscillator drive. Un- 
expectedly, the rim velocity initially undershoots, recov- 
ers from the undershoot exponentially with a character- 
istic time constant and eventually decays logarithmically 
in the long time limit. The time constant tends to diverge 
as temperature is decreased below 15 mK. 
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